RAL-93-092 
November 1993 



Improved QCD sum rule estimates of the higher twist contributions 
to polarised and unpolarised nucleon structure functions 

o 

Q 

| Graham G. RossQ 

Department of Theoretical Physics, University of Oxford, 
; Oxford, UK. 

CO 
CN 

CNJ \ and 

m 

Q\ ■ R.G. Roberts 



a,: 



X 



Rutherford Appleton Laboratory, 



q; Chilton, Didcot 0X11 OQX, UK. 

43 



Abstract 

We re-examine the estimates of the higher twist contributions to the integral of g\, 
the polarised structure function of the nucleon, based on QCD sum rules. By including 
corrections both to the perturbative contribution and to the low energy contribution we 
find that the matrix elements of the relevant operators are more stable to variations 
of the Borel parameter M 2 , allowing for a meaningful estimate of the matrix elements. 
We find that these matrix elements are typically twice as large as previous estimates. 
However, inserting these new estimates into the recently corrected expressions for the 
first moments of g\ leads to corrections too small to affect the phenomenological analysis. 
For the unpolarised case the higher twist corrections to the GLS and Bjorken sum rules 
are substantial and bring the estimate of Aqc d from the former into good agreement with 
that obtained from the Q 2 dependence of deep inelastic data. 
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The recent measurements of the polarised nucleon structure functions g[ ,n,d [0 , @ , ||] raise 
the question of the consistency of the three measurements Ip^a = Jo dxgi ,n,d (x). In particular 
since the three measurements are at different values of < Q 2 >, namely 10.7, 2 and 4.6 GeV 2 
for p, n and d, it is crucial to understand the Q 2 dependence of the first moment of g\ in order 
to test the Bjorken sum rule |4j or to extract an estimate of the nucleon's spin content, Aq, 
consistent with all three experiments. 

The higher order corrections to the leading twist expressions have been calculated to 0(a 2 ) 
for the non-singlet quantities and to 0(a s ) for the singlet contribution but it is the power 
corrections from the higher twist operators which have recently ||, || been shown may play 
an important role in a consistent analysis of the data. The magnitudes of the reduced matrix 
elements of the relevant higher twist operators U s , U NS , V s , V NS were extracted from a QCD 
sum rule calculation by Balitsky, Braun and Kolesnichenko (BBK) and used in the analysis 
of ref 0. The aim of the present paper is to sharpen the results of BBK by re-examining the 
computation of <C JJ S ' NS <C V S,NS including a contribution to the perturbative QCD 
side of the sum rule which was dropped in the Borel transformation and explicitly retaining 
the continuum term to the nucleon pole on the low-energy side of the sum rule. This leads to a 
significant improvement in the stability of the extracted value of the reduced matrix elements 
and hence to a more reliable estimate of the higher twist contribution to the integral of g±. 
Moreover the estimated values of the matrix elements are significantly larger than previously 
but since the coefficient of the twist three piece in the first moment moment of g\ has very 
recently been corrected ||], it turns out that the net higher twist contribution to the moment 
is minimal. The improvement in stability and increased magnitude of the matrix elements is 
also true for the unpolarised case and we find that the correction to the Gross-Llewellyn Smith 
(GLS) sum rule|| is sufficient to affect the extracted value of Ajjg substantially. 

Following the procedure of BBK we consider the quantity if we are interested in the 

operators U S ' NS , given by 

r» = i 2 f dx e*» f dy {T[n(x)Up(y)fj{0))) W 
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where rj is the nucleon current. Expanding in powers of -y gives 



A V " ln a (-£j) + B ln(^) + C(— 2 ) ln(-^) + C{\) + £>(i 
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The coefficients A, ..,D may be read off from eq(8) of BBK corresponding to the QCD evaluation 
of T^, including non-perturbative effects due to QCD condensates. The next step is to Borel 
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transform the coefficient of 2^75 in eq(@) which gives 
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(3) 

which differs from eq(lf ) of BBK since we have explicitly carried out the p 2 integration of the 
p-ln(^5-) contribution instead of neglecting itQ. 

To complete the sum rule the quantity T^fj?) must also be determined in terms of the nucleon 
and continuum contributions. Balitsky et al.fTj use the form 
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(4) 



(p 2 — m 2 ) 2 (p 2 — m 2 ) 

where the first term is the pure nucleon pole contribution and the single pole term is added to 
allow for the interference of the pole term with a continuum contribution. The Borel transform 
of the coefficient of 2^75 in eqn([|) is 
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and the QCD sum rule results from equating eq(§) to eq(|D to give 
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In this paper we are particularly concerned to estimate of the errors in determining the 
operator matrix elements from the QCD sum rules. Thus we will consider in detail the effects 
of each of the terms in this expansion and the inclusion of further terms in the Borel expansion 
in (M 2 ) n /n\. On the rhs of eq(|2]) the first term not included is ~ 1/p 6 which, after Borel 
transformation contributes a term ~ 1/2!M 2 to the rhs of eq(||). The coefficient of this term 
should be no bigger than the coefficient of Y if the perturbative expansion in l/(n\(M) n ) is 
acceptably convergent. By adding such a term and fitting its coefficient we may check this over 
some range of M 2 and hence establish over what range of M 2 (if any) the sum rule converges. 

The corrections to the lhs of eq(H) are somewhat more difficult to determine. Further 
resonance contributions with mass m 2 R > m 2 can be included by adding a resonance pole 

IGeV 2 , this term vanishes. However the 



2 Balitsky et al argue that with the natural choice, fj, 2 
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Borel transform integrates over p 2 so this is, at best, an approximation. Since the effect of the term can be 
explicitly included we choose not to make this approximation. 
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(m R — p 2 )^ 1 - After Borel transformation this gives an additional term M 2 e ( - m2 ~ m R^ M ' 2 Y on 
the lhs of eq(^]). However the dominant correction to the lhs of eq(^) is not expected to come 
from a nucleon resonance excitation but from the (N + nir) continuum which has a threshold 



at E — ym 2 + m 2 , very close to the nucleon pole. As far as we know this contribution has not 
been considered explicitly even though it is potentially very large. However, we will demonstrate 
that this contribution does not substantially degrade the accuracy with which we can determine 
the operator matrix elements provided the resonance term discussed above is added to eq(|]). 
The reason is that the ttN contribution is well described by the terms proportional to X in 
ec L@ together with the resonance term proportional to Y. To demonstrate this we will estimate 
the contribution due to the Ntc intermediate state, fig. 2. This gives the term 

T 91nnK f /(g 2 ) 2 [~2(p + qW + « U >] A m 
* N m 2 J (2vr) 4 (g 2 -m 2 )(m 2 -(p + g) 2 ) 2 U 

Here we have included a form factor f(q 2 ) needed to describe the rjNir coupling far from 
the pion mass shell (on shell we take (77 | Ntt) = X p g n NN)- The result is insensitive to the 
particular choice of f(q 2 ) provided it provides convergence for large q 2 ; here we choose f{q 2 ) = 
q 2 /((q 2 ~ m2 r)(l ~~ q 2 /0.7GeV 2 )) i.e. chosen to vanish as q 2 — > and to be the same as the 
nucleon electromagnetic form factor for large spacelike q 2 . (Since f{q 2 ) is needed for spacelike 
q 2 we have chosen a form factor with no singularities in the spacelike region.) After taking 
the Borel transform of the coefficient of p^frys and multiplying by M 2 e m2//M2 /2A 2 we find the 
resulting contribution to the lhs of eq@ is accurately reproduced by the form 

a + f3M 2 + 7 M 2 e (m2 - m * )/M2 (8) 

with a ~ 0.03 < U >, (5 ~ 0.15 < U >, 7 ~ 0.2 < U >, and m\ w 2m 2 . The origin of this 
form is easy to understand. The ttN singularity starts at (m n + m) in the energy plane and, 
after Borel transformation, gives a contribution proportional to the integral of e (m ~ m ' over 
m, with an appropriate weighting factor. To a good approximation this may be approximated 
by the sum of two exponentials e ( m2 - m i)/ M2 anc [ e (m 2 -m 2 R )/M 2 w }th mi ~ ( m?r -f ra). Taylor 
expanding the first term in leads to eq(|8|). Thus we see that the ttN contribution may, to a 
good approximation, be included in the continuum and interference terms and gives a relatively 
small correction to the determination of the operator matrix elements. Using similar methods 
suggests that the contribution to constant term in eq@ from N + mc intermediate states will 
also be quite small. 



3 



Let us now turn to the phenomenological analysis and first consider the analysis of ref@. 
We have four operators to consider U s , U NS , V s , V NS the latter two having p^7s replaced by 
Si/, a ^n,uPaPfj,luj5 where A, S stand for symmetric and antisymmetric combination of indices. 
So we have four QCD sum rules with coefficients A — > E given by BBK eqs(8,9). In extracting 
a value of <C U ^> from the sum rule, BBK retained only the first two terms on the lhs and also 
set C—0. They obtained <C U ^> at each value of M 2 by applying (1 — M 2 -^^) to the rhs. The 
stability of this estimate of <C U 3> relies upon the rhs of eq(|6]) being approximately linear over 
a range of M 2 . The solid lines in fig.l are the BBK values for the matrix elements extracted by 
this procedure and the variation with M 2 reflects the importance of higher derivatives in M 2 , 
making it difficult to arrive at a reliable estimates for <C U, V ^>. 

Let us see how the situation changes with the addition of the proposed continuum term 
M 2 e (m 2 ~ m |)/M 2 Y and the term proposal to C in eq(Bl). As we will demonstrate these are 
necessary to describe adequately the M 2 behaviour of the rhs over a reasonable range. We 
first fit the rhs of eq(^) by the form on the lhs with four parameters C [/ >, X, F and 
m 2 R (the last corresponding to an effective threshold for the continuum) over a wide range 0.5 
GeV 2 < M 2 < 2.0 GeV 2 . Then, with m 2 R fixed (typically around 2 GeV 2 ), we perform a three 
parameter fit over each interval of 0.2 GeV 2 in M 2 . The dashed line in fig. 2 shows the resulting 
values of -C U ^> and we immediately see a marked improvement in the stability compared 
to neglecting the continuum term which results in a more meaningful estimate of the matrix 
elements. At the same time we notice that the actual value obtained is considerably different. 
At M 2 = 1 GeV 2 for example the new estimate is roughly twice the old value. If instead we 
choose to describe the continuum by e ( m2 - m2 R )/ M2 y (j. e . a double resonance pole) there is a 
similar improvement in stability and the extracted values of the matrix elements are within 
15% of those in fig.l. 

Using these results we will now try to sharpen the estimates of the higher twist contributions 
further by considering possible additional terms to the rhs of eq(|6|). We will consider the 
sensitivity of <C U, V ^> to an additional term ~ ^ on the rhs of eq(|]), i.e. a term ~ ^ifi 
on the rhs of eq(g). Thus we may have an additional term 6 2Af2 on the lhs of eqn (6) which 
can be used in fitting the rhs and help determine the range of M 2 over which the sum rule is 
valid. As a result of fitting the curves (with /z 2 = 1 GeV 2 ) in fig.l with the constraint that the 
coefficient of the additional term is no bigger than the coefficient Y of the continuum term lead 
to an acceptable range 0.8 < M 2 < 1.75 GeV 2 . 
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In fact, the magnitude of the continuum is correlated, as one might expect, with the value 
of the hadronic continuum parameter s on the rhs. Our results correspond to Sq = 2.25 GeV 2 
as in refs[[7|, [10], 11 1. The dependence on sq is weak however and careful fitting reveals that 
this dependence is absorbed by the explicit continuum contributions, leaving the magnitude of 
<C U, V ^> practically invariant as Sq varies in the range 1.8 to 5 GeV 2 . 

For the above range in M 2 , the resulting uncertainties in the values of <C U, V ^> from the 
fitting procedure are comparable to the uncertainties expected from varying fi 2 in the range 
0.33 to 3 GeV 2 . The values of the reduced matrix elements obtained are 

< U s >= 0.046 ± O.OlOGeV 2 < U NS >= 0.317 ± O.OWGeV 2 
m 2 < V s >= -0.292 ± OMOGeV 2 m 2 « V NS >= 0.605 ± 0.030GW 2 (9) 

From the values in eq(|9|) we compute the coefficients a p and a n of the 1/Q 2 contributions 
to the integrals of g\ and g™, using the corrected formulas of ref[§ 



a p + a n - 
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< U s > -\m 2 < V s > 
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(10) 



which gives 

a„ = -0.029 ±0.002, a„ = -0.002 ± 0.002. (11 
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The errors in eq([5])are based simply on the small variation of the values of <C U, V ^> with M 
and fi 2 and are typically 5%. These errors are in addition to the underlying uncertainties arising 



from the factorisation assumption for the vacuum condensates which are typically 20% [12 
Thus a realistic estimate of the errors in eq(|Tl~D is more like 0.010. 
The integrals I p , n ,d can be written 

I p = h + I 8 + I + a p /Q 2 
In = -h + h + h + a n /Q 2 



I 8 + I + (a p + a n )/Q 2 (12) 



where 
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Using the measured values of the polarisation asymmetries from refs |T], H, HI, the vam es 
of I p ,n,d a ^ values of Q 2 = 10.7, 2, 4.6 GeV 2 were extracted in ref |5| and determined to be 
0.134±0.012, -0.023±0.005, 0.041±0.016 respectively. The estimates for the coefficients a p , a n 
from our improved QCD sum rule analysis, eq(|i~l~D, when inserted into eqs (p!2]JT3| ) yield estimates 
for the nucleon spin content (for F/D = 0.575 ± 0.016) 

Aq = 0.24 ±0.11 fromp 
= 0.53 ±0.07 fromn 

= 0.27 ±0.15 fromd (14) 

The value of Aq obtained from the relatively low Q 2 neutron data from SLAC is still out of 
line with the values obtained form CERN on the proton and deuterium. If the neutron higher 
twist coefficient a n had come out large and positive, around 0.04 or so, then Aq would decrease 
a value close to the proton and deuteron estimates. In fact the analysis of Ellis and Karliner|| 
used such a value based on the BBK[|7| estimates of <C U, V ^> but with the (now known to 
be) incorrect formulas for twist three contribution to the first moment. Thus despite the fact 
that we claim considerably larger estimates for the matrix elements, the corrected formulas|| 
lead to a small neutron correction, thereby ruling out higher twists as a way of reconciling the 



three experiments. Bag model estimates for a n give a zero value JT3[ . 

The improvement to our understanding of the QCD sum rule estimates of -C U, V ^> has 
led to a more meaningful determination of the higher twist corrections to the integrals of the 
polarised structure function g\. We recall that two corrections to the sum rule — the evaluation 
of the p-ln(^j) contribution on the rhs and the inclusion of the continuum contribution — 
resulted in more stable estimates of the reduced matrix elements. It is therefore natural to 
ask if similar corrections apply to other sum rules e.g. the matrix elements <C O s,NS ^> which 
determine the 1/Q 2 corrections to the GLS|§ and Bjorken (Fi) ||14|| sum rules flO| . The first 



correction does not apply since hi(p-) terms cancel for the unpolarised operator p^p 1 but the 
continuum correction should be included. 

Fig. 3 shows the corresponding determinations of <C O s,NS In fig. 3(a) we note the strong 
variation of s ^> with M 2 , again indicating the inadequacy of the nucleon pole terms alone 
in eq(|J). Carrying out an analogous fitting procedure as in the polarised case, the resulting 
estimate for s ^> is far less sensitive to the value of M 2 indicated by the dashed line in 
fig.3(a). Interestingly, the estimated magnitude of both matrix elements <C O s ' NS 3> increases 
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when a realistic continuum term is included. In particular, the estimate used in the analysis of 
Chyla and Kataev[I5| was that of ref[]IU 



< s >= 0.33±0.16GW 2 (15) 

which led to a value of Aj^g extracted from the data on xF 3 (x, Q 2 = 3 GeV 2 on the GLS sum 
rule of 

A^L = 318 ± 23(stat) ± 99(syst) ± 62(twist)MeV. (16) 
If Iqls is the measured value of the GLS sum rule, then 

Igls _ a s (Q 2 ) 8 « s » 
3 7T 27 Q 2 1 } 

and we see that an increased estimate of <C s 3> leads to a lower value of Ajjg. We estimate 
the larger and more precise value 

< s >= 0.53 ± OMGeV 2 (18) 

which leads to 

Ag_ = 232 ± 23(stat) ± 99(syst) ± 17(twist)MeV r (19) 

which is more in accord with estimates got from studying the Q 2 dependence of deep inelastic 
data@, Ag^ = 230 ± 55 MeV. A gain, the error in eq(|18|) does not include the intrinsic 
uncertainty (~ 20%) associated with the factorisation assumption; including this raises the 
final error in eq(^) from 17 to 45 MeV. 

In summary, we have shown that there are corrections to QCD sum rules which have not been 
included in previous determinations of the relevant reduced matrix elements. These corrections 
lead to significant improvement in the stability of the extracted values with respect to the 
range of the Borel parameter M 2 . As an example, we have applied these corrections to the 
case of the polarised structure function gi of the nucleon and found that the size of the the 
higher twist contributions is now determined with better precision. We have considered in 
some detail the corrections to the sum rule. The Nn correction, potentially very significant, 
has been shown to have little effect on the determination of the reduced matrix elements. The 
remaining contributions are under control for a reasonable range of M 2 allowing for a realistic 
determination of the error. Our final estimates for the matrix elements, both for polarised 
and unpolarised structure functions are more reliable and, moreover, significantly larger than 
previous estimates. Nevertheless because of the recent corrections to the term multiplying the 



twist three contribution, the resulting impact on the phenomenology of moments of g^ ,n,d is 

MS 



reduced. For the analysis of the GLS sum rule however, the resulting value of is more in 



line with other deep inelastic scattering phenomenology. 
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Figure Captions 

Fig. 1 Values (in GeV 2 ) of the matrix elements <C U, V 3> derived from fitting the rhs of eq(P) 
in bins of 0.2 GeV 2 in M 2 with just the two terms on the lhs (solid lines) and together 
with the continuum term M 2 e ( - m2 ~ m R^ M2 V(dashed lines). 

Fig. 2 it N contribution, eq(|7]), to the correlation function. 

Fig. 3 Values (in GeV 2 ) of the matrix elements <C O s ' NS ^ derived from fitting the rhs of eq(39) 
of BK|10[ with the nucleon double and single pole terms only (solid lines) and together 
with a continuum term (dashed lines), as in fig. 1. 
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